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Chapter  1 


INTRODUCTION 

Standard  Base  Supply  System 

The  Standard  Base  Supply  System  (SBSS)  is  an  auto¬ 
mated  Inventory  accounting  and  control  system  designed  to 
provide  timely  support  to  base  level  activities.  The  system 
uses  the  UNIVAC  1050-11  computer  for  storage  and  maintenance 
of  records  and  for  the  generation  of  management  reports. 

The  system  as  a  whole  consists  of  both  manual  components  and 
interfacing  computer  programs.  The  system  consists  of  four 
major  functional  processes --item  accounting,  accounting  and 
finance,  file  maintenance,  and  management  reporting  (24 i 1-2). 
When  using  the  UNIVAC  1050-11  on-line  computer  system,  all 
records  affected  by  a  given  transaction  are  updated  at  the 
time  of  transaction  input.  The  system  provides  remote 
devices  in  work  centers  on  base  to  insure  timely  communica¬ 
tions  between  man  and  machine  (24il-5). 

The  SBSS  was  designed  to  help  speed  the  flow  of 
materials  from  many  soxirces  to  the  base  user  organizations. 
The  SBSS  inventory  is  stocked  by  material  from  Air  Force 
depots,  the  General  Services  Administration  (GSA) ,  Defense 
Supply  Agency  (DSA)  and  through  local  procurement  (19il). 

The  SBSS  is  driven  by  user  demand  actions. 


1 


The  S3SS  automatically  provides  outputs  to  other 
management  systems.  Those  outputs,  such  as  requisitions, 
go  to  vrtiolesale  suppliers  such  as  Air  Force  Logistics  Com¬ 
mand  (AFLC)  and  the  Defense  Logistics  Agency  (DLA) .  The 
S5SS  imposes,  Air  Force  wide,  standard  organization,  pro¬ 
grams  and  procedures.  The  organization  is  designed  for  ease 
of  customer  support  and  efficiency  of  internal  operation 
(23il-6),  The  entire  SBSS  is  an  extension  of  the  basic 
supply  requirements  to  order,  receive,  store,  and  issue 
property. 


Inventory 

Inventories  are  used  throughout  the  world  in  many 
organizations,  both  civilian  and  military.  Inventories  can 
be  classified  in  many  ways.  The  following  classification 
is  one  convenient  wayi 

1.  Production  inventoiriesi  raw  materials,  parts 
and  components  vAiich  enter  the  orgsmizatlon's  product  or 
service  in  the  production/ tramsformatlon  processes. 

These  may  be  either  special  items  manufactured  to  user 
specifications  or  standard  items  purchased  "off  the 
shelf." 

2.  MRO  inventories!  maintenance,  repair  and  oper¬ 
ating  supplies  consumed  in  the  production/ transformation 
processes  but  which  do  not  become  part  of  the  product  or 
service, 

3.  In-process  inventories!  semifinished  products 
found  at  various  stages  in  the  production  operation. 

4.  Finished -goods  inventories!  completed  products 
ready  for  shipment  and/or  use  £l4!l89J. 

Regardless  of  their  classification,  inventories 
serve  the  primary  purpose  of  decoupling  successive  stages 
in  the  production-distribution-consumption  chain  (I3!l). 
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Inventories  allow  production  decisions  and/or  supplier 
demand  decisions  to  be  made  independent  of  supplier  procvire* 
ment  decisions  where  demand  is  Independent  or  near  independ¬ 
ent. 

Inventory  Control 

Inventory  control  is  a  vital  element  in  the  manage¬ 
ment  of  system  resources.  Development  of  analytical  tech¬ 
niques  and  computer  capability  have  combined  to  transform 
Inventory  control  into  a  critical  function  requiring  profes¬ 
sional  managerial  skills  (14tl87). 

The  significance  of  the  Inventory  control  function 
Is  clearly  demonstrated  by  the  level  of  resources  Involved, 
Studies  have  sho%m  that  a  commercial  firm’s  Inventory  com¬ 
monly  constitutes  anywhere  from  15  to  25  percent  of  Its  total 
Invested  capital  (14il87),  In  Its  everyday  operations,  In¬ 
ventory  control  Is  often  a  computerized  operation  within  a 
carefully  defined  and  controlled  structural  framework 
(14 1 22 5).  Sound  management  In  the  area  of  Inventory  control 
Is  vital.  Management  must  try  to  balance  the  various  risks 
of  low  Inventory  with  Its  associated  risk  of  stock  out, 
production  halts,  back  orders i  and  high  Inventory  with  Its 
risk  of  high  carrying  cost  and  Increased  risk  of  obsoles¬ 
cence  (23 1142).  To  achieve  this  memagement  objective  a 
forecast  of  future  demand  must  be  made.  This  forecast  Is 
based  on  uncertain  economic  cycles,  customer  demands,  and 
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advances  In  technology.  In  light  of  this  uncertain  environ¬ 
ment,  inventory  control  must  answer  how  much  to  order  and 
when  to  order  it  (l3i2). 

Generally  speaking,  three  basic  systems  have  been 
employed  in  controlling  Inventories i  (1)  the  cyclical  order¬ 
ing  system,  (2)  fixed  order  quantity  system,  and  (3)  mate¬ 
rial  requirements  planning  system  (14 i 194).  The  first  sys¬ 
tem  is  a  time-based  system  involving  scheduled  periodic 
reviews  of  the  stock  level  of  all  Inventoiy  items.  The  sec¬ 
ond  type--the  fixed  order  quantity  system--is  based  on  the 
order  quantity  factor  rather  than  on  a  time  factor.  The 
material  requirements  planning  system  concept  provides  a 
way  of  looking  at  the  management  of  production  inventories 
in  a  dependent  demand  envixponment  (14il94),  Regardless  of 
the  system  used,  management  must  have  an  accurate  assessment 
of  what  will  happen  in  the  future  in  order  to  make  its  deci¬ 
sions  (13 1 2). 

Forecasting 

In  almost  any  activity,  having  perfect  knowledge  of 
the  future  would  be  an  immense  help.  Rather  than  the  cer¬ 
tainty  of  perfect  knowledge,  the  inventory  manager  finds 
himself  in  an  uncertain  environment.  He  can  do  little  to 
change  that  environment  so  he  seeks  management  techniques 
which  will  aid  in  reducing  the  levels  of  uncertainty.  All 
forecasting  techniques  have  in  common  the  basic  goal  of 
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reducing  the  uncertainty  in  one's  expectation  of  the  future 
(4il).  As  Robert  L.  Sims  says,  "It  should  be  obvious  that  I 

to  forecast  is  to  be  wrong.  Only  in  the  most  fortuitous 
cases,  ax\d  these  seldom  occur,  can  a  forecaster  guess  the 
future  exactly  C2lil]."  Today,  the  economics  of  spiraling 
prices  and  austere  budgets  make  reliable  estimates,  or  fore¬ 
casts,  especially  critical  (3i22).  Forecasting  for  and 
acquiring  spare  parts  is  an  Important  facet  of  the  material 
support  system.  Forecasts  are  used  to  develop  dependable 

I 

systems  and  to  formulate  decisions  concerning  the  techno-  j 

logical  and  economic  feasibility  of  retaining  a  system  in  ' 

! 

use  (111  15).  The  exponential  acceleration  of  spares  cost 
provides  a  compelling  reason  for  basing  requirements  corapu- 

! 

ration  upon  the  best  available  forecasting  methodology  i 

(3i22). 

Standard  Base  SuodIv  Forecasting  System 

The  Air  Force  Logistics  Command  (AFLC),  throughout 
its  Air  Logistics  Centers  (ALCs)  or  depots,  buys  and  dis¬ 
tributes  centrally  procured  items  of  supply  used  on  mis¬ 
siles,  aircraft,  and  other  equipment.  Each  depot  buys  and 
stocks  items  in  specific  federal  supply  classes  (FSCs)  and 
serves  as  the  primary  source  of  supply  for  Air  Force  bases 
and  other  activities  (24 • 11-1). 

The  SBSS  is  driven  by  user  demand  actions.  When  a 
base  needs  an  item  a  requisition  is  submitted  to  SBSS.  If 
the  item  needed  is  not  currently  in  stock,  a  request  for 
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that  item  is  placed  on  order.  The  request  may  be  filled 
through  the  normal  operating  cycle.  In  some  instances,  spe¬ 
cial  orders  may  be  needed.  Orders  for  resupply  of  the  SBSS 
are  logged  as  due-ins  euid  are  maintained  until  the  material 
is  received  ( 19 i 1) , 

Repair  parts  and  supplies  with  expendability, 
repairability,  recoverability  category  (ERRC)  designator 
XB3  will  be  stocked,  using  a  variable  economic  order  quan¬ 
tity  (EOQ)  stockage  concept.  An  XB3  item  is  an  expendable 
item,  t.e.  an  item  consumed  in  use  or  incorporated  into 
another  assembly.  The  EOQ  model  is  applicable  when  the 
quantity  of  items  ordered  arrives  in  the  inventory  at  one 
point  in  time  and  when  the  demand  for  the  item  has  a  con¬ 
stant,  or  nearly  constant,  ratei  and  the  cost  of  the  item 
is  relatively  low  (2i497),  The  SBSS  requirements  determina¬ 
tion  may  be  divided  into  the  range  model  and  a  depth  model. 
The  range  model  refers  to  the  procedure  used  to  determine 
if  an  item  is  to  be  stocked  at  the  base  level.  The  depth 
model  is  used  if  the  item  is  to  be  stocked  to  determine  how 
much  to  order  and  when  to  order  it.  The  depth  model  is  based 
on  the  classic  EOQ  formula t 

EOQ  -  y  2DA/IP  (1.1) 

where,  D  *■  annual  demand  rate 

A  ■  cost  per  order  (currently  used  figure,  $4,54) 

I  "  annual  Inventory  carrying  rate  (currently  used 
figure,  26  percent) 

P  ■  item  unit  price  (19 i2) 


This  formula  balances  the  cost  of  ordering  with  the  cost  of 
holding  inventory.  The  derived  order  quantity  will  minimize 
these  total  variable  costs  (I9i2). 

The  Air  Force  application  of  the  EOQ  concept  pro¬ 
duces  a  variable  requisition  objective  by  determining  the 
EOQ  in  consideration  of  the  number  of  demands,  daily  demand 
rates,  and  a  stockage  priority  code.  The  purpose  of  the 
variable  EOQ  stockage  concept  is  to  prevent  the  premature 
stockage  of  items  based  on  erratic  or  indefinite  demand  pat¬ 
terns  (24 j 11-2),  The  basic  EOQ  model  has  been  enhanced  to 
better  meet  the  veuriable  stockage  objectives  (VSO)  of  the 
SBSS« 

EOQvso  =  V  2«DDR*VS0«A/(I.P)  (1.2) 

where,  A  =  cost  per  order  ($4.54) 

I  =  annual  inventory  carrying  rate  (26  percent) 

P  =  item  unit  price 
DDR  =  daily  demand  rate 

^  VSO  =  number  of  days  in  the  EOQ  computation  (19i3) 
Another  important  concept  incorporated  into  the 
SBSS  is  that  base  supply  levels  are  managed  on  a  reorder 
point  basis.  That  is,  each  time  the  stock  position  of  an 
item  reaches  or  falls  below  the  established  reorder  level, 
some  type  of  supply  action  must  be  taken  to  bring  stocks 
up  to  quantities  authorized  (24 t 11-2).  The  reorder  point 
is  a  combination  of  the  order  and  shipping  time  quantity 
(OSTQ)  and  the  safety  level  quantity  (SLQ) .  The  OSTQ  is 
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that  quantity  required  by  a  base  to  permit  uninterrupted 
replacement  from  the  external  supply  source.  The  OSTQ  is 
given  by I 


OSTQ  =  DDR*OST 


(1.3) 


where,  DDR  =  daily  demand  rate 

OST  =  average  order  and  shiptime  in  days,  based  on 
the  item  source  and  priority  (19i3) 

The  safety  level  quantity  is  a  variable  quantity,  in  units, 

computed  to  provide  protection  from  stockouts  during  the 

reordering  process.  The  SLQ  is  given  by« 


SLQ  =  C  V  3. OSTQ  (1.4) 

where,  C  =  the  safety  factor  (typically  set  to  1,  which 
implies  an  84  percent  service  effectiveness, 
given  the  assumption  of  normally  distributed 
demands) 

3  =  lead  time  demand  variance/mean  ratio 
OSTQ  =  order  and  shiptirae  quantity  (19i3) 

The  reorder  point  (RP)  is  given  byi 


RP  =  OSTQ  +  SLQ  (19i3)  (1.5) 


The  maximum  desired  inventory  position  is  called  the  requi¬ 
sition  objective  (RO) i 


RO  =  INT  (E0Qy3q  +  OSTQ  +  SLQ  +  0.999)  (1.6) 


vrfiere,  INT  =  integer  result  of  the  following  computation 
EOQvso  “  enhanced  EOQ  model 

OSTQ  “  order  and  shiptime  quantity 
SLQ  ■  safety  level  quantity 
0.999  “  factor  added  to  the  EOO  level  to  adjust  to 
the  next  highest  unit  C19i4) 
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It  can  be  seen  that  each  of  the  three  variables  in  the 
requisition  objective  computation  are  a  function  of  the 
daily  demand  rate.  The  daily  demand  rate  is  the  forecast 
measurement  for  the  SBSS.  If  inaccurate  estimates  for  DDR 
are  made,  errors  may  result  (19i3). 

'.Vhen  an  item  is  ordered  the  daily  demand  rate  is 
revised.  The  SBSS  maintains  status  on  the  cumulative  recur¬ 
ring  demand  (CRD)  and  the  date  of  first  demand  (DOFD) .  Each 
time  an  item  is  ordered,  the  number  of  units  ordered  is 
added  to  the  CRD.  The  revised  DDR  is  them 

DDR  =  CRD/max  (180,  Current  Date  -  DOFD)  (1.7) 

The  minimum  of  180  days  usage  is  assumed  so  as  not  to  over¬ 
stock  items  that  have  been  recently  added  to  the  stockage 
list  (19i3).  The  number  of  demands  for  each  ERRC  code  XB3 
item  will  be  recorded  in  six-month  increments  up  to  a  total 
of  three  increments  or  18  months*  i.e.,  the  current  and  two 

s. 

past  six-month  increments.  When  18  months*  demands  have 
been  accumulated,  the  oldest  six-month  increment  will  be 
dropped  and  a  new  six  months'  accumulation  will  begin*  The 
criteria  for  establishing  demand  levels  for  EOQ  items  (ERRC 
code  XB3)  will  vary  dependent  upon  the  stockage  priority 
code,  but  will  not  be  less  than  three  demands  per  year 
(24tll-3). 
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Problem  Statement 


The  literature  search  has  shown  that  a  significant 
error  variance  exists  in  the  forecasting  method  employed  in 
the  Standard  Base  Supply  System  \&ien  compared  to  the  actual 
d«nand  experienced.  The  current  forecasting  system  used  by 
the  SBSS  is  a  combination  of  exponential  smoothing  and  mov¬ 
ing  averages.  The  literature  suggests  that  this  technique 
may  not  provide  the  best  estimate  of  future  demand. 

A  critical  measurement  used  in  SBSS  is  forecasted 
demand.  A  forecasting  approach  »dilch  reduces  forecast  error 
has  potential  impact  beyond  its  effect  upon  economic  order 
quantity  (EOQ)  determination.  Accurate  demand  data  in  the 
context  of  the  SBSS  affect  average  Inventory  levels  and 
holding  costs,  and  improve  order  and  shipping  time  (lead- 
time),  thereby  impacting  safety  stock  levels  and  premium 
transportation  requirements  for  priority  back  orders. 

This  thesis  explores  whether  a  multiple  model  fore¬ 
casting  technique  can  Improve  demand  forecasts. 

Research  Oblectives 

The  objectives  of  this  thesis  arei 

1,  Evaluate  the  multiple  model  forecasting  tech¬ 
nique  as  an  alternative  to  the  Standard  Base  Supply  System 
forecasting  technique. 

2.  Determine  which  forecasting  approach  will  pro¬ 
vide  the  most  accurate  estimate  of  future  usage. 


3.  Determine  the  effectiveness  of  implementing  the 
proposed  system. 

4,  Recommend  actions  based  upon  the  results 
obtained . 

Research  Justification 

This  thesis  applies  the  research  of  Mitchell  and 
Garland  at  the  wholesale  level  of  demand  to  retail  level 
forecasting.  Those  authors  used  simulated  and  limited  actual 
data  to  conclude  that  the  multiple  model  forecasting  tech¬ 
niques  actually  forecasted  more  accurately  than  the  single 
method  currently  used  in  the  Air  Force's  wholesale  level 
D062  system  (I3i42).  This  thesis  had  the  sponsorship  of  the 
Air  Force  Logistics  Management  Center  (AFLMC) ,  Gunter  Air 
Force  Station,  in  a  continued  effort  to  Improve  present 
Inventory  control  methods. 

Scope 

This  thesis  will  utilize  a  reseaurch  data  base  main¬ 
tained  by  AFLMC  and  containing  two  and  a  half  yeaurs  of  data 
from  Dover  Air  Force  Base,  Delaware,  The  data  base  contains 
approximately  15,000  line  items  (National  Stock  Numbers). 

The  data  were  tested  using  several  different  forecasting 
techniques.  The  results  were  compared  to  results  obtained 
by  previous  research  on  the  SBSS. 
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Chapter  2 


LITERATURE  REVIE’.^ 


Background 

Major  Richard  Lombardi  of  the  Stocking  Policies 
Division,  Air  Force  Logistics  Management  Center,  stimmarized 
the  current  dilemma  of  logisticians  faced  with  forecasting 
requirements  to  meet  mission  demands i 

The  Standard  Base  Supply  System  (SBSS)  demand  fore¬ 
casting  is  an  vinorthodox  prediction  scheme  which  may 
lead  to  suboptimal  forecasting.  This  is  attributed  to 
the  manner  in  vrtiich  the  daily  demand  rate  (DDR)  is  up¬ 
dated  by  six  month  adjustments  in  the  SBSS  ctimulated 
recurring  demand  (CRD)  and  the  date  of  first  demand 
(DOFD)  fields.  The  net  effect  of  those  adjustments  is 
to  convert  the  forecasting  to  a  form  of  exponential 
smoothing  in  which  the  value  of  the  smoothing  constant 
will  vary  depending  upon  date  of  demands.  This  varia¬ 
tion  may,  depending  upon  when  demands  are  experienced, 
provide  more  weight  to  past  demands  at  the  expense  of 
current  demands,  thus  dampening  the  influence  of  cur¬ 
rent  requirements 


Multiple  Model  Forecasting  Method 

In  order  to  choose  an  effective  forecasting  system, 
the  forecaster  must  decide  upon  a  model  vrtiich  is  most  appro¬ 
priate  given  the  conditions  that  exist  at  the  time  of  the 
forecast.  A  variety  of  forecasting  models  are  available 
from  which  the  forecaster  c^ul  choose.  The  selection  of  a 
forecasting  method  depends  on  many  factors i  the  context  of 
the  forecast,  i.e.,  the  natiore  of  the  decision  environment 
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ac  the  time  of  the  forecast,  the  relevance  cind  availability 
of  historical  data,  the  desired  accuracy  of  the  forecast, 
and  the  time  period  to  be  covered  (13 ill).  Due  to  these 
factors,  certain  types  of  forecasting  techniques  are  better 
suited  to  a  particular  demand  pattern  than  others.  Since 
in  practice  we  are  seldom  faced  with  a  ptire  "classical" 
demand  pattern,  the  forecaster  should  consider  several 
methods  for  forecasting  a  single  demand  pattern  before 
selecting  the  model  best  suited  to  the  actual  demand  pat¬ 
tern. 

The  multiple  model  forecasting  method  used  in  this 
thesis  research  incorporates  forecasting  techniques  recom¬ 
mended  by  a  variety  of  experts.  These  techniques  are  moving 
average,  single  exponential  smoothing,  double  exponential 
smoothing,  and  adaptive  exponential  smoothing.  The  computer 
was  programmed  with  these  simple  forecasting  strategies. 
Through  the  process  of  simulation,  the  computer  will  select 
the  one  best  strategj  to  forecast  an  item  at  a  given  moment 
in  time.  Whichever  strategy  best  projected  the  most  recent 
completed  period  is  the  one  the  computer  uses  to  project 
the  future  (22 i 3). 

Focus  Forecasting 

A  new  approach  to  Inventory  control  using  the  con¬ 
cept  of  declining  computer  costs  per  computation  is  advo¬ 
cated  by  Bernard  T.  Smith  in  his  book.  Focus  Forecasting! 
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Computer  Techniques  for  Inventory  Control  (22).  Focus  fore¬ 
casting  as  implemented  by  Smith  used  a  series  of  simple 
forecasting  approaches.  It  uses  the  powerful  computation 
simulation  ability  of  the  computer  to  pick  the  one  forecast¬ 
ing  strategy  that  will  work  best  for  the  one  item  for  the 
next  period  (22ixii).  Focus  forecasting  goes  through  all 
the  computations  every  time  it  forecasts  a  demand. 

Several  parallel  research  efforts  have  been  accom¬ 
plished  using  various  forecasting  techniques  at  the  wholesale 
level  on  the  AFLC  D062  inventory  control  system  for  expend¬ 
able  items  and  the  AFLC  D041  reparable  asset  management  sys¬ 
tem.  A  summary  of  these  theses  follows. 

Garland  and  Mitchell 

The  purpose  of  their  study  was  to  determine  if  a 
multiple  model  forecasting  technique  could  forecast  demand 
more  accurately  than  the  model  currently  used  in  the  Air 
Force  Logistics  Command  D062  system  for  expendable  items. 
Simulated  and  actual  data  were  used  to  check  the  results. 

The  methods  utilized  in  the  multiple  model  technique  were 
an  eight-quarter  moving  average*  a  four-quarter  moving  aver¬ 
age,  exponential  smoothing,  adaptive  smoothing,  a  least 
squares  fit  and  a  ratio  of  change  between  years  method. 
Results  were  compared  in  terms  of  mean  absolute  deviation 
adjusted  to  show  percentage  change  in  accuracy  compared  to 
the  D062.  The  statistical  test  used  for  comparison  was  the 
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t-test  for  matched  pairs.  This  test  indicated  approximately 
a  l7  percent  improvement  in  accuracy  using  either  simulated 
or  real  historical  data. 

Christensen  and  Schroeder 

This  research  effort  compared  the  D041  Single  Moving 
Average  forecasting  method  used  to  forecast  reparable  gen¬ 
erations  of  recoverable  items  with  the  Box  and  Jenkins*  time 
series  analysis  forecasting  methods.  Five  artificially  gen¬ 
erated  stochastic  processes  were  used  to  model  the  possible 
reparable  generations  observed  in  practicei  (1)  a  Poisson 
process  with  a  constant  mean,  (2)  a  Poisson  process  with  a 
decreasing  mean,  (3)  a  Poisson  process  with  an  alternating 
linear  mean,  and  (4)  a  process  whose  values  are  the  sine 
function  of  the  output  of  a  Poisson  process.  The  research 
concluded  that  the  D041  forecasting  method  made  unbiased 
forecasts  for  the  Poisson  process  with  a  constant  mean  and 
the  sine  function,  but  made  biased  forecasts  for  the  other 
three  processes.  Time  series  analysis  forecasting  methods 
were  only  used  to  make  forecasts  for  the  processes  that 
were  found  to  be  biased  using  the  D041  forecasting  method. 
Time  series  analysis  forecasting  methods  made  unbiased  fore¬ 
casts  for  the  processes  whose  means  were  linearly  increasing, 
linearly  decreasing,  and  alternating  linearly. 
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Brantley  and  Loreman 

This  research  effort  examined  the  forecasting  tech¬ 
nique  used  in  the  Air  Force  Logistics  Command  D041  reparable 
asset  management  system.  It  was  hypothesized  that  the  mean 
of  the  absolute  values  of  the  D041  forecast  error  in  prac¬ 
tice  was  equal  to  zero.  This  hypothesis  could  not  be 
rejected.  It  was  further  hypothesized  that  another  time 
series  forecasting  technique  (exponential  smoothing)  would 
also  yield  an  error  distribution  with  a  mean  of  zero.  This 
hypothesis  could  not  be  rejected  eitherj  moreover,  the  vari¬ 
ances  for  the  exponential  smoothing  forecast  error  distribu¬ 
tions  were  less  than  the  D041  forecast  error  distributions 
for  all  four  lead  times  examined  (one,  two,  three  and  four 
quarters) . 


Fischer  and  Gibson 

This  thesis  effort  examined  the  application  of 
exponential  smoothing  to  forecast  demand  for  economic  order 
quantity  items.  The  research  team  stated  that  manag«nent  of 
EOQ  items  required  the  use  of  a  demeuid  forecasting  technique 
to  estimate  future  demand  for  the  purpose  of  establishing 
stock  levels.  The  study  compared  the  effectiveness  of  four 
forecasting  models  that  could  be  used  at  base  level.  The 
moving  averages  method  and  single,  double  and  triple  expo¬ 
nential  smoothing  were  evaluated  using  22  months  of  demand 
history  for  a  random  sample  of  34  EOQ  Items  stocked  at  a 


base  consolidated  supply  activity.  Four  statistical  error 
measures  were  used  to  compare  the  accuracy  of  the  forecasts 
generated  by  the  models  for  the  items.  The  study  concludes 
that  the  methods  were  not  significantly  different  in  fore¬ 
casting  EOQ  items  in  terms  of  accuracy,  stability  and  bias. 
The  study  team  recommended  a  reporting  system  be  considered 
-hat  would  permit  the  study  of  the  application  simultaneously 
of  various  forecasting  techniques  to  the  management  of  an 
item  and  the  factors  affecting  item  demand  patterns, 

Patterson 

The  purpose  of  this  research  was  to  look  into  alter¬ 
nate  approaches  to  forecasting  demand  for  expendable  items 
in  the  SBSS.  The  forecasting  models  studied  were  single, 
double  and  adaptive  exponential  smoothing.  Sample  items 
were  selected  from  the  actual  demand  records  from  Dover  AFB, 
An  analysis  of  the  various  models  was  performed  by  means  of 
a  computer  program  written  for  each  model.  The  forecasting 
models  were  compared  on  the  basis  of  forecast  error  as  meas¬ 
ured  by  the  mean  absolute  deviation  (MAD) ,  The  forecast 
error  for  the  model  currently  used  by  the  SBSS  was  also 
measured.  The  research  study  concluded  that  the  single 
exponential  smoothing  method,  using  small  smoothing  con¬ 
stants,  produced  the  lowest  error  rate. 
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Conclusion 

The  literature  search  has  revealed  that  the  current 
forecasting  system  used  by  the  SBSS  may  not  provide  the 
best  estimate  of  futiire  demand.  Additionally,  the  literature 
search  has  shown  the  SBSS  forecasting  method  exhibits  a  sig¬ 
nificant  error  variance  when  compared  to  actual  demand. 

These  facts  lead  to  several  research  questions  which  this 
research  effort  will  explore. 

Research  Questions 

1.  Is  there  a  significant  error  variance  in  the 
multiple  model  forecasting  technique  when  compared  to  pre¬ 
vious  research? 

2.  '.Vhlch  of  the  forecasting  techniques  selected 
for  this  study  exhibits  the  smallest  error  variance,  i.e., 
v^ich  technique  most  accurately  forecasts  demand? 

Hypothesis 

The  multiple  model  forecasting  technique  performs 
as  well  as,  or  better  than,  single  model  techniques  to  which 
it  is  compared. 
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Chapter  3 


METHODOLOGY 


Introduction 

The  methodology  used  in  this  research  effort 
requires  data  to  be  gathered  and  then  £malyzed  using  multi¬ 
ple  model  forecasting  techniques.  The  results  obtained  from 
the  multiple  model  forecasting  technique  will  be  compared 
to  the  results  of  techniques  used  In  previous  research. 

Using  this  procedure  recommendations  on  the  performance  of 
the  multiple  model  forecasting  technique  can  be  made. 

Model  Selection 

The  current  Standard  Base  Supply  System  (SBSS)  Is  a 
mlxtiire  of  moving  averages  and  exponential  smoothing*  It 
Is  a  moving  average  since  the  Dally  Demand  Rate  (DDR)  Is  a 
quotient  of  cumulative  demands  over  time  divided  by  the 
length  of  time  for  which  these  demands  are  accumulated  (see 
equation  1.7).  It  Includes  exponential  smoothing  by  the 
way  updates  to  Cumulative  Recurring  Demand  (CRD)  and  Date 
of  First  Demand  (DOFD)  are  performed  (19 i 6).  It  Is  not 
apparent,  however,  that  the  current  SBSS  approach  provides 
the  best  estimate  of  future  demand.  An  approach  to  fore¬ 
casting  which  reduces  forecast  error  should  provide  a  more 
accurate  estimate  of  future  demand.  The  techniques 
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considered  in  this  thesis  are  moving  average,  single  expo¬ 
nential  smoothing,  double  exponential  smoothing,  adaptive 
exponential  smoothing,  and  multiple  model  forecasting. 

Moving  Average 

\’ihen  demand  for  an  item  does  not  have  a  rapid  growth 
or  seasonal  characteristics,  a  moving  average  can  be  useful 
for  forecasting.  The  moving  average  model  would  be  expected 
to  perform  well  in  normal  constant  demand  situations.  This 
method  asstames  that  the  data  generating  process  constitutes 
a  time  series i 


=  D  +  e 


(3.1) 


where,  F^.  *  forecasted  demand  for  period  t 

D  =  average  demand  over  time 
e  *  random  error  variable  with  a  mean  of  zero 
as  constant  variance  over  time  (20 i 103) 

The  moving  average  technique  may  be  described  asi 


F^l  -  °t  Vl  Vz  Pt-ntl  (3.2) 

N 


t-rr*-! 

1  y  r 

«  1-t 


where , 


t-rrM 


i-t 


N 


forecast  for  the  next  period,  t+l 

actual  demand  at  time  t,  t-1,  t-2  ...  t-ntl 

the  number  of  observations  used  in  the 
coverage 
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The  effect  of  moving  average  forecasting  depends  on 
the  number  of  observations  used  in  the  coverage  (N) .  If  N 
is  large,  the  equal  weighting  given  to  each  period  is  small, 
and  random  fluctuations  have  little  effect  on  the  foirecast. 

If  N  is  small,  the  model  is  more  sensitive  to  changes  in 
demand  (I2tl7). 

Single  Exponential  Smoothing 

In  using  exponential  smoothing,  three  pieces  of  data 
are  needed i  the  most  recent  forecast,  the  actual  demand 
that  occurred  for  that  forecasted  period,  and  a  smoothing 
constant  (Ot).  The  equation  for  a  single  exponential  smooth¬ 
ing  forecast  isi 

♦  Q((D^  -  Fj)  (3.3) 

where,  ■  the  exponentially  smoothed  forecast  for 

periodVl 

■  the  exponentially  smoothed  forecast  for  the 
prior  period  t 

■  the  actual  demand  in  the  prior  period  t 

0(  "  the  desired  response  rate,  or  smoothing 

constant  (I3i5d) 

This  equation  states  that  the  new  forecast  is  equal  to  the 
old  forecast  plus  an  adjustment  proportional  to  the  differ¬ 
ence  between  the  previous  forecast  and  the  actual  experience. 
The  closer  CX  is  to  1,  the  more  the  new  forecast  will  incor¬ 
porate  an  adjustment  for  the  error  in  the  upcoming  forecast. 
The  closer  Of  is  to  0,  the  less  sensitive  the  new  forecast 
will  be  to  the  error  in  the  prior  forecast  (I3i58). 
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Double  Exponential  Smoothing 


Double  exponential  smooching  is  an  extension  of 
single  exponential  smoothing*  It  is  usually  applied  to 
items  that  exhibit  a  trend  pattern.  Both  single  and  double 
smooched  values  lag  actual  data  when  a  trend  exists,  the 
difference  between  single  and  double  smoothed  values  can  be 
added  to  the  single  smooched  values  and  thereby  adjust  for 
trend  (19 ill).  The  single  exponential  model,  equation  (3.3), 
can  be  rewritten  asi 


St  “  Sc.i  +  a(D^  -  (3.4) 

%rtiere,  “  exponentially  smooched  dem£md  in  period  i. 

The  apparent  trend  component  of  the  error  factor  of  our 
forecast  isi 


I 


-  S 


t-l 


(3.5) 


This  trend  component  can  also  be  smoothed  or  averaged  over 
time  usings 

*  ■'t-i  ♦ 

The  double  exponentially  smoothed  forecast  or  the  smoothed 
average  forecast  including  the  trend  component  is  the  simple 
smooched  average  plus  the  smootl^  trend  component  as  cor¬ 
rected  for  lag  in  T^  by  the  term  (l-0')/C^i 


S  +  — ^T 
t  at 


(3.7) 
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In  order  to  use  this  approach  only  three  data  values 
and  a  smoothing  constant  are  required.  Use  of  this  model 
should  be  based  on  assumptions  about  the  demand  pattern. 

If  the  time  series  has  a  trending  average  demand  rate,  the 
double  exponential  smoothing  model  produces  an  accurate 
estimate  of  demand  (I2il9). 


Adaptive  Exponential  Smoothing 

The  adaptive  smoothing  techniques  adjust  Of  over 
time  based  on  the  size  of  the  forecasting  error.  The  adap¬ 
tive  model  utilized  here  is  the  same  as  in  single  exponen¬ 
tial  smoothing  with  one  exception.  The  value  for  Oi  is 
derived  from  the  equationi 


a 


t+i 


M. 


(3.8) 


where, 


D  -  F 
t  t  t 

P  ■  smoothing  constant 
and  I  I  denotes  absolute  values  (16t54). 

The  characteristic  of  not  specifying  a  value  for  Oi 
is  attractive  when  thousands  of  items  require  forecasting. 
Also,  this  method  can  change  the  value  of  Cf  when  changes 
in  the  pattern  of  the  data  have  made  the  initial  Oi  value 
no  longer  appropriate  (16i53}. 
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Multiple  Model  Forecasting 

The  model  that  most  accurately  forecasts  demand  will 
be  selected  from  all  models  considered.  In  order  to  do 
this,  each  demand  history  will  be  broken  Into  three  periods i 
a  base  (or  start-up)  period,  a  test  period,  emd  a  prediction 
period. 

The  base  (or  start-up)  period  Is  the  historical  data 
base  required  for  each  model.  The  length  of  the  base  period 
varies  from  model  to  model. 

After  the  forecasting  start  point  Is  established, 
each  model  forecasts  demand  for  a  given  test  period.  The 
forecasts  are  compared  to  actual  demand  for  the  test  period. 
The  model  with  the  smallest  variation  from  the  actual  demand 
is  selected  to  forecast  demand  for  the  next  prediction 
period. 

This  is  a  recurring  process  in  that  the  forecasting 
horizon  defines  a  new  prediction  period)  the  past  prediction 
period  becomes  the  current  test  period)  the  past  test  period 
Joins  the  base  period  data  base)  and  older  data  are  dropped 
from  the  base  period  (13il2). 

The  multiple  model  method  employed  in  this  research 
is  shown  in  Appendix  D,  subroutine  "Multi."  The  actual 
demand  for  the  four  previous  forecasting  periods  is  added 
and  stored  as  a  new  variable.  The  values  forecasted  for 
these  same  four  periods  by  each  technique  are  added  and  this 
total  is  subtracted  from  the  total  for  the  actual  demand. 
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This  process  Is  repeated  for  each  forecasting  technique. 

The  technique  which  displays  the  lowest  variance  from  the 
actual  demand  Is  selected  by  the  multiple  model  technique 
as  the  technique  It  will  use  to  forecast  demand  for  the 
next  period. 

Accuracy  of  Forecasting  Models 

Demand  for  an  Item  Is  generated  through  the  inter* 
action  of  a  number  of  factors.  This  interaction  is  extremely 
complex.  Due  to  this  complexity*  all  forecasts  will  contain 
some  error. 

In  many  forecasting  situations  accuracy  Is  treated 
as  the  main  criterion  for  selecting  a  forecasting  method. 

In  spite  of  this  fact,  little  systematic  work  has  been  done 
to  develop  a  framework  for  measuring  and  evaluating  accuracy- 
related  Issues.  One  of  the  difficulties  In  dealing  with  the 
criterion  of  accuracy  In  forecasting  situations  Is  the 
absence  of  a  single  universally  accepted  measure  of  accuracy 
(16iS69).  When  discussing  forecast  errors  It  Is  appropriate 
to  distinguish  between  sources  of  error  and  the  measurement 
of  error. 

Sources  of  error  can  be  classified  as  either  bias 
or  random.  Bias  errors  occur  when  a  consistent  "mistake” 

Is  made.  Random  errors  can  be  defined  as  those  that  cannot 
be  explained  by  the  forecast  model  employed. 
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Several  of  the  common  terms  used  to  measure  the 
degree  of  error  are  mean  absolute  percentage  error  and  mean 
absolute  deviation.  Additionally,  measures  of  bias  may  be 
used  to  indicate  the  amount  of  any  positive  or  negative  bias 
in  the  forecast  (9i24l), 


Mean  Absolute  Percentciae  Error 

One  approach  used  to  detect  forecast  error  is  to 
take  the  absolute  value  of  each  individual  percentage  error 
to  obtain  the  mean  absolute  percentage  error  (MAPE) . 

Percentage  error  is  a  relative  meastire  of  accuracy 
and  is  defined  asi 


PE^  - 


-  F. 


(100) 


(3.9) 


where,  PE^  “  percentage  error  for  any  time  period  t 
“  actual  value  for  time  period  t 
■  forecast  value  for  time  period  t  (I6i570) 

The  equation  for  computing  the  MAPE  over  N  time 


periods  is  given  byi 


MAPE 


N 


(3.10) 


where,  MAPE  ■  mean  absolute  percentage  error 

PEj^  ■  percentage  error  for  a  time  period  (I5i570) 

However,  the  MAPE  itself  does  not  always  give  a  good 

basis  of  comparison  as  to  the  gains  in  accuracy  made  by 
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applying  a  specific  forecasting  method,  MAPE  indicates  how 
well  the  forecasting  model  fits  the  data.  In  this  way,  it 
aids  in  evaluating  how  accurate  the  model  is  (16i568). 

Mean  Absolute  Deviation 

The  most  straightforwaird  measure  of  error  is  mean 
absolute  deviation  (MAD) .  It  is  computed  using  the  differ¬ 
ences  between  the  actual  demand  and  the  forecasted  demand 
without  regard  to  sign.  Since  it  is  the  mean  deviation,  the 
sum  of  the  absolute  deviations  is  divided  by  the  number  of 
data  points.  Stated  in  equation  formi 

t  |“t  -  '■t 

MAD  »  -Sil -  (3.11) 

N 

where,  D^  •  actual  demand  for  period  t 

■  forecasted  demand  for  period  t 
N  *  total  number  of  periods 

and  I  I  denotes  the  absolute  value  (9i242). 

The  mean  absolute  deviation  ignores  whether  the  forecast  is 

greater  or  less  than  actual  demand  2uid  simply  measures  the 

magnitude  of  difference. 

Bias 

Bias  indicates  the  directional  tendency  of  forecast 
errors.  If  the  forecasting  technique  used  consistently  over¬ 
estimates  actual  demand,  bias  will  be  positlvei  consistently 
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cause  the  bias  to  be 


underestimating  actual  demands  will 
negative.  Bias  is  defined  asi 

n 

I -  “t)i 

aa  9  -  A*1 _ _ 


Where,  »  forecasted  demand  for  period  t 
Dj.  *•  actual  demand  for  period  t 
N  »  total  number  of  periods  (li33i) 

Conclusion 

All  three  techniques  described  here  are  used  to 
evaluate  multiple  model  forecasting.  The  measures  of  error 
In  the  forecasting  techniques  are  used  to  determine  which 
forecasting  method  will  provide  the  most  accurate  estimate 
of  demand.  By  using  the  forecasting  technique  which  is  the 
most  accurate,  the  error  variance  in  the  Standard  Base  Supply 
System  can  hopefully  be  reduced. 


► 
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Chapter  4 


DATA  COLLECTION  AND  ANALYSIS 


Introduction 

The  purpose  of  this  chapter  is  to  identify  the 
source  of  data  used  in  the  analysis  and  describe  the  selec¬ 
tion  procedure  used  to  obtain  a  saunple  of  279  National 
Stock  Numbers  (NSNs).  Additionally,  this  chapter  will  pre¬ 
sent  an  analysis  of  that  data. 

Data  Collection 

As  mentioned  in  Chapter  1,  this  research  had  the 
sponsorship  of  the  Air  Force  Logistics  Management  Center 
(AFLMC),  Gunter  AFS,  Alabama.  Two  computer  tapes  were  pro¬ 
vided  which  contained  Stand2u:d  Base  Supply  System  data  from 
Dover  AFB,  Delaware,  for  a  two  and  one-half  year  period 
from  October  1976  to  March  1979,  One  tape  contained  item 
record  data  for  33  stock  classes,  ranging  from  the  1005  to 
the  53l5  stock  class.  The  total  number  of  stock  classes  was 
115  containing  a  total  of  10,607  individual  National  Stock 
Numbers.  The  remainl.ng  tape  contained  transaction  history 
data  for  these  same  NSNs.  A  Fortran  prograun  (Appendix  A) 
was  developed  udiich  aggregated  these  records  by  stock  class} 
total  items  in  each  class i  and  total  recurring  demands. 
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The  next  step  was  selecting  the  stock  classes  to  be 
used  In  the  actual  analysis.  A  previous  study  by  Patterson 
(19)  testing  alternative  forecasting  techniques  for  the 
Standard  Base  Supply  System  (SBSS)  used  stock  classes  28,  47, 
59,  and  66.  It  was  the  researchers*  Intention  to  select  a 
sampling  of  these  same  classes  In  order  to  make  a  comparative 
analysis  of  results.  However,  since  the  data  only  Included 
through  the  53  class,  two  additional  alternative  stock 
classes  were  chosen.  A  Fortran  program  (Appendix  B)  was 
used  to  prioritize  the  data  into  groups  of  high  activity 
(frequency  of  demand)  and  quantity  demanded.  With  this 
Information,  the  authors  were  able  to  Identify  the  51  and 
53  stock  classes  as  having  characteristics  that  would  lend 
themselves  to  alternative  forecasting  methods  and,  partic¬ 
ularly,  to  the  multiple  model  method.  Table  I  shows  a  cap- 
sullzed  summary  of  the  activity  and  demands  for  the  four 
stock  classes  selected. 


Table  1 

Profile  of  Preliminary  Data 


Stock  Class 

Total  NSN's 

Total 

Recurring  Demands 

28 

386 

18,196 

47 

1313 

84,939 

51 

2407 

126,571 

53 

2331 

842,654 

30 


Due  to  the  large  frequency  and  quantity  of  demand 
data  and  the  limited  computer  resources,  a  sampling  of  this 
population  was  taken.  Recurring  demands  were  tabulated  from 
the  transaction  history  tape  (see  Fortran  program  in 
Appendix  C)  and  aggregated  into  130  weeks  or  data  points  for 
each  selected  item. 

Prior  to  the  selection  of  the  final  sample,  the  data 
were  scanned  for  any  obvious  outliers.  As  in  Patterson's 
research,  one  stock  number  (47  stock  class)  had  6,160  units 
requested  in  one  period  and  zero  demands  for  the  remaining 
129  periods.  This  item  was  dropped  from  the  sample.  A  28 
stock  class  item  was  also  discarded  because  data  were  fovind 
to  be  available  only  for  the  first  one  and  one-half  years  of 
the  survey  period. 

Table  2  is  a  summary  of  the  final  stock  classes 
used  in  testing  the  alternative  forecasting  methods. 

Table  2 

Profile  of  Stock  Classes  Selected  for  Analysis 


Stock  Class 

Total  NSN's 

Total 

Recurring  Demands 

28 

99 

8,137 

47 

99 

14,131 

51 

31 

21.659 

53 

50 

234,689 
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In  sviinmary«  the  data  used  for  the  actual  analysis 
consisted  of  centrally  procured  items  identified  by  the 
Expendabilityi  Repairability*  Recoverability  Category  (ERRC) 
of  XB3.  A  total  sample  of  279  line  items  was  selected 
based  on  high  activity  and  a  relatively  high  frequency  of 
reorder  periods.  An  inherent  bias  in  the  data  collection 
was  the  researchers'  attempt  to  identify  those  stock  numbers 
which  exhibited  characteristics  which  would  accommodate  the 
time  series  forecasting  techniques  to  be  exercised. 

This  research  addressed  the  applicability  of  multi¬ 
ple  model  fozrecasting  as  an  alternative  to  the  standard  base 
supply  forecasting  system.  In  order  to  test  this  hypothesis, 
the  data  collected  by  the  methods  described  were  processed 
by  use  of  a  Fortran  IV  computer  program  (see  Appendix  0). 

Computer  Output 

For  each  stock  number  in  each  stock  class,  the 
actual  units  demanded,  the  \inits  forecasted  by  each  fore¬ 
casting  technique  used,  and  the  units  forecasted  by  the  mul¬ 
tiple  model  forecasting  technique  were  shown  on  the  output 
(see  Appendix  E) .  The  mean  absolute  deviation  (MAD)  com¬ 
paring  the  forecasted  demand  and  the  actual  demand  was  com¬ 
puted  (see  equation  3.11).  Additionally,  a  bias  (equation 
3.12)  and  mean  absolute  percentage  error  (MAPS,  equation 
3.10)  were  computed.  These  values  were  necessary  to  complete 
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the  statistical  tests  required  to  determine  the  accuracy  of 
the  various  forecasting  methods. 


Data  Base  Comparison 

In  order  to  attain  the  objective  of  evaluating  the 
multiple  model  forecasting  method.  It  Is  necessary  to  show 
that  a  similar  data  base  exists  with  which  to  compare  our 
results.  In  his  research,  Patterson  obtained  a  mean  abso¬ 
lute  deviation  for  single  and  double  exponential  smoothing, 
using  an  alpha  of  .2,  among  others  (19i9).  A  comparison  of 
the  means  for  these  two  forecasting  techniques  will  test  the 
data  base  used  to  see  If  It  Is  essentially  the  same  as  that 
used  by  Patterson  In  his  research.  The  precise  stock  numbers 
used  by  Patterson  were  unknown  to  us. 

To  test  this  assumption,  95  percent  confidence  Inter¬ 
vals  were  computed.  It  Is  necessary  to  calculate  a  value  for 
the  standard  deviation  for  Patterson's  research,  since  none 
was  given.  When  forecast  errors  are  normally  distributed, 
the  MAD  Is  related  to  the  standard  deviation  byi 


1  Standard  Deviation 


where,  TT  ■  3,1416 
or  conversely. 


(4,1) 


1  MAD  ■  0.8  Standard  Deviation  (9i242)  (4.2) 
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The  central  limit  theorem  tells  us  that  for  almost  all  popu¬ 
lations  the  sampling  distribution  of  the  sample  mean  Is 
approximately  normal  when  the  sample  size  Is  sufficiently 
large  (I8i202).  Since  the  sample  size  used  by  Patterson  Is 
large  (200),  normality  can  be  assumed,  and  a  value  can  be 
calculated  for  the  standard  deviation. 

Using  this  value  for  the  standard  deviation  of  both 
samples,  a  95  percent  confidence  Interval  was  coastructed. 

To  construct  a  confidence  Interval,  an  unbiased  estimator  of 
the  difference  In  the  two  sample  means  Is  required.  This 
estimator  Is  denoted  by  Di 

D  -  Y  -  X  (4.3) 


where,  Y  ■  mean  of  observations  from  sample  2 

Y  ■  mean  of  observations  from  sample  1  (I8t3l2) 

Additionally,  an  unbiased  estimator  for  the  standard  devia¬ 
tion  of  D  Is  required I 


S^(D)  -  S^(Y)  S^(X) 


(4.4) 


trtiere,  S^(Y)  ■ 


n 


S^(X)  -  (I8i3.6) 


The  fozmiula  for  calculating  the  two-sided  confidence  Inter¬ 
val  Is  I 


L  -  5  -  a(l-0(/2)S(D) 

U  -  C  »(l-a/2)S(lJ)  (18i316) 


(4.5) 
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The  results  of  the  confidence  interval  calculations 


are  shown  In  Table  3: 


Table  3 

95%  Confidence  Interval  for  Data  Base  Comparison 


Model 

Lower 

Limit 

Mean 

MAD 

Upper 

Limit 

Stock  Class  28 

Single  Exponential 

.6282 

2.11 

2.5918 

Double  Exponential 

.8670 

2.36 

2.8931 

Stock  Class  47 

Single  Exponential 

.5062 

2.06 

2.2938 

Double  Exponential 

.2405 

2.47 

2.6396 

Since  the  means  fell  within  the  limits*  we  have  concluded 
that  for  practical  purposes  our  data  base  and  that  used  by 
Patterson  were  essentially  the  same. 

Data  Analysis 

The  results  of  the  computer  analysis  of  the  data 
are  shown  in  Table  4.  In  three  of  the  four  stock  classes, 
on  the  average,  single  exponential  smoothing  consistently 
displayed  the  lowest  mean  absolute  deviation.  This  is  con¬ 
trary  to  what  was  expected.  The  data  were  selected  from 
the  various  stock  classes  based  on  large  amounts  of  demand 
occurring  in  the  most  time  periods.  Essentially  the  "deck 
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Table  4 


Results  of  Computer  Analysis 


Mean  MAD 

Mean  Bias 

Mean  MAFE 

Stock  Class  28 

Single  Exponential 

2.1115 

• 

.6691 

88.7535 

Double  Exponential 

2.3659 

• 

.3165 

114,7806 

Adaptive  Exponential 

2.9618 

.5948 

178.4009 

4 -Term 

2.3750 

• 

.2396 

117.1407 

8 -Terra 

2.2967 

• 

.3491 

111.1047 

Multiple  Model 

2.7200 

- 

.2081 

152.5846 

Stock  Class  47 

Single  Exponential 

2.0682 

• 

1.0852 

213.0647 

Double  Exponential 

2.4756 

- 

.1933 

278.1608 

Adaptive  Exponential 

3.8648 

1.3457 

388.1105 

4-Terra 

2.5769 

• 

.2905 

265.9312 

8-Term 

2.3863 

• 

.4338 

251.1996 

Multiple  Model 

2.8990 

.6206 

328.0684 

Stock  Class  51 

Single  Exponential 

13.6679 

m 

1.6520 

615.5107 

Double  Exponential 

14.8354 

.8196 

781.0071 

Adaptive  Exponential 

18.1501 

3.8259 

1081.6932 

4 -Terra 

14.4877 

m 

.3390 

710.0943 

8-Tertn 

13.9992 

m 

.6066 

677.3004 

Multiple  Model 

16.1558 

'1.8196 

894.6500 

Stock  Class  53 

Single  Exponential 

98.2697 

• 

3,9447 

4023.2286 

Double  Exponential 

107.5422 

14.2225 

5036.5975 

Adaptive  Exponential 

114.3334 

15.5944 

5373.4977 

4 -Terra 

103.0422 

1.0954 

4245.1873 

8 -Term 

97.7790 

m 

2.0539 

4073.0743 

Multiple  Model 

107.2436 

8.0751 

4762.1244 

36 
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was  stacked"  in  favor  of  the  time  series  techniques.  The 
multiple  model  technique  did  not  produce  the  most  accurate 
estimate  of  demand,  nor  did  it  even  tie  with  the  most  accu¬ 
rate  method.  In  the  case  of  the  51  and  53  stock  classes, 
the  MAD  for  all  techniques  considered  rises  considerably. 

This  indicates  the  problems  encountered  in  using  time  series 
analysis  on  data  that  are  highly  erratic  and  sparse  (many 
zero  entires).  The  MATE  reflects  this  highly  erratic 
nature.  The  relatively  large  values  for  the  MAPEs  indicate 
abrupt  changes  in  the  demands  in  the  data  ( see  Appendix  £) . 
There  are  numerous  periods  in  %diich  zero  demand  is  incurred 
and  suddenly  a  large  change  in  demand  occurs  and  abruptly 
drops  off  to  zero  again. 

Conclusion 

Through  the  use  of  confidence  intervals  it  was  shown 
that  the  data  base  used  in  this  research  was  essentially  the 
sane  as  that  used  by  Patterson  in  his  previous  research.  The 
Fortran  IV  program  produced  results  that  indicated  single 
exponential,  and  not  multiple  model  forecasting,  was  the 
most  accurate  of  the  techniques  considered.  The  implica¬ 
tions  and  conclusions  wliich  have  been  drawn  from  this  will 
be  further  discussed  in  the  next  chapter. 
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Chapter  5 


CONCLUSIONS  AND  RECOMMENDATIONS 


Introduction 

This  research  was  directed  at  exploring  the  hypoth¬ 
esis  i  the  multiple  model  forecasting  technique  performs  as 
well  as,  or  better  than,  single  model  techniques  to  which 
it  is  compared.  The  research  addressing  this  h3rpothesis  led 
to  several  conclusions,  vdiich  can  be  associated  with  the 
four  objectives  of  this  thesis. 

Conclusions 

The  first  objective  of  this  research  was  to  evaluate 
the  multiple  model  forecasting  technique  as  an  alternative 
to  the  Standard  Base  Supply  System  forecasting  technique. 
Closely  related  to  this  was  the  second  objective,  to  deter¬ 
mine  which  forecasting  approach  will  provide  the  most  accu¬ 
rate  estimate  of  future  demand. 

The  multiple  model  forecasting  technique  and  the 
additional  techniques  were  processed  by  the  use  of  the  For¬ 
tran  IV  program  referred  to  In  Chapter  4.  Table  4  in  Chap¬ 
ter  4  reflects  the  results  of  the  data  analysis.  As  can  be 
seen,  multiple  model  forecasting  did  not  produce  the  most 
accurate  estimate  of  demand.  Based  on  the  mean  MAD  computed 
after  the  computer  analysis  was  completed,  multiple  model 
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consistently  ranked  next  to  last  In  accuracy.  On  three  out 
of  the  four  stock  classes,  single  exponential  smoothing  was 
ranked  as  the  most  accurate.  On  the  fourth  class,  stock 
class  53,  single  exponential  smoothing  was  the  second  most 
accurate  method,  following  closely  behind  8>term  moving 
average. 

The  bias  is  an  indicator  of  the  directional  tendency 
of  the  forecast  errors  for  each  forecasting  technique.  In  a 
perfect  forecasting  model  the  bias  would  be  expected  to  be 
zero,  i.e.,  the  random  errors  due  to  overforecasting  would 
be  offset  by  random  errors  due  to  underforecasting.  Again, 
referring  to  Table  4,  in  three  out  of  the  four  stock  classes 
analyzed,  multiple  model  forecasting  consistently  on  the 
average  overforecasts  the  demand.  The  implications  of 
overforecasting  demand  come  into  play  when  the  holding  cost 
of  an  item  is  considered.  Holding  cost  is  important  to  any 
manager  concerned  with  any  form  of  inventory.  Holding  cost, 
along  with  order  costs  and  demand  rates,  makes  up  the  eco¬ 
nomic  order  quantity  (EOQ),  The  multiple  model  forecasting 
technique  used  in  this  research  would  cause  the  holding  cost 
to  rise,  if  only  slightly,  due  to  its  consistent  over¬ 
estimating  of  demand  and  subsequent  over- stocking  of  inven¬ 
tory.  Single  exponential  smoothing,  as  used  in  this 
research,  consistently  under-estimates  demand  and  therefore 
again  becomes  a  more  acceptable  forecasting  technique  ^en 
compared  to  others,  including  multiple  model. 


Finally,  the  mean  absolute  percentage  error  (MAPE) 
of  the  multiple  model  technique  consistently  ranks  as  the 
next  to  least  accurate.  The  MAPE  indicates  how  well  the 
technique  in  question  fits  the  data.  Therefore,  the  multi¬ 
ple  model  forecasting  used  in  this  research  does  not  fit 
the  data  as  well  as  other  methods.  Single  exponential 
smoothing  fits  the  data  best,  l.e..  It  consistently  e^diib- 
Ited  the  lowest  MAPE,  The  high  MAPEs  encountered  in  this 
research  indicate  the  highly  erratic  nature  of  the  SBSS. 

An  example  of  the  erratic  nature  of  the  demand  data  is  shown 
in  Appendix  E.  Additionally,  it  can  be  seen  that,  as  a  gen¬ 
eral  rule,  trend  la  not  apparent  In  SBSS  Items.  Therefore, 
methods  of  forecasting  that  did  adjust  for  trend  did  not 
reduce  the  forecast  error.  When  an  abrupt  change  occurred 
in  the  data,  both  double  and  adaptive  exponential  smoothing 
reacted  to  the  change  as  if  it  were  a  trend.  However,  the 
data  usually  abruptly  returned  to  zero  and  a  trend  did  not 
develop.  As  a  result,  these  forecasting  methods  were  not 
among  the  roost  accurate  of  the  methods  tested. 

The  third  objective  of  this  research  was  to  deter¬ 
mine  the  effectiveness  of  Implementing  the  proposed  system. 
In  order  to  demonstrate  a  significant  cost  savings  to  the 
Air  Force,  there  must  be  a  significant  Increase  in  fore¬ 
casting  accuracy.  The  results  Indicate  that  the  multiple 
model  forecasting  technique  used  here  would  not  support 
this  significant  Increase  In  forecasting  acciiracy. 


Therefore,  any  cost  savings  attributable  to  implementing  the 
multiple  model  forecasting  technique  are  in  doubt.  It  also 
becomes  apparent  that  the  time  series  forecasting  models 
used  in  this  research  did  not  yield  any  significant  improve¬ 
ment  in  the  forecasting  of  demand  for  SBSS  items. 

The  final  objective  of  this  research  was  to  recom¬ 
mend  actions  baaed  upon  the  results  obtained.  The  following 
section  will  address  this  objective. 

Recoramendati bns 

Implied  throughout  this  study  is  the  fact  that 
Inventory  models  are  used  to  assist  managers  in  determining 
when  items  should  be  ordered  and  how  large  the  order  should 
be  so  as  to  minimize  variable  costs  and  insure  that  adequate 
stockage  is  on-hand  to  meet  mission  requirements.  Using  the 
EOQ  model  as  a  depth  model,  however,  has  two  inherent 
assumptions,  i.e.,  the  demand  activity  remains  relatively 
continuous  and  the  items  are  relatively  low  in  cost.  Refer¬ 
ring  to  Table  5,  we  note  that  this  model  was  intended  for 
use  under  the  conditions  which  are  present  in  Quadrant  I. 

Table  5 

Inventory  Characteristics 


Inventory  Item  Cost 
HI  LO 


HI 

II 

I 

Item 

Demand 

LO 

III 

IV 
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Quadrant  1  exhibits  the  ideal  conditions  for  EOQ,  i.e.,  a 
high  frequency  of  demand  with  continuous  demcuid  behavior. 
Scarcity  of  resources  and  other  fiscal  environmental  factors 
have  made  the  conditions  within  Quadrant  II  not  at  all  un* 
common.  Average  item  costs  have  been  increasing  signifi¬ 
cantly  in  recent  years.  The  £0Q  model,  while  not  optimal, 
can  handle  these  conditions  fairly  well  as  long  as  the 
demand  patterns  remain  relatively  continuous.  In  addition 
to  costs,  technology  factors  have  increased  reliability  on 
many  items  resulting  in  lower  demand  frequency.  This  has 
the  effect  of  moving  many  items  into  Quadrant  III  and  IV. 

An  example  would  be  vacuum  tubes  which  are  now  replaced  by 
solid  state  technology  resulting  in  a  higher  Mean  Time  Be¬ 
tween  Failure,  Quadrant  III  and  IV  contain  characteristics 
of  low  demand  for  both  low  and  high  value  items.  The  math¬ 
ematical  assumptions  inherent  within  the  EOQ  model  were  not 
meant  to  address  the  conditions  of  low  demand,  especially 
the  sparse  Ixjmpy  conditions  found  frequently  in  the  SBSS 
data.  Examples  of  this  demand  data  are  found  in  Appendix  E. 

We  conclude  that  the  EOQ  model  is  an  inappropriate 
depth  model  for  items  characterized  by  erratic  and  sparse 
demand  patterns.  Our  recommendation  is  to  divert  further 
efforts  from  improving  SBSS  forecasting  techniques  to  direct¬ 
ing  research  efforts  in  a  search  for  alternative  decision 
models  idiich  will  be  more  appropriate  for  the  SBSS  depth  of 
stockage  decision. 
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APPENDICES 
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APPENDIX  A 

DATA  RETRIEVAL  PROGRAM 
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ion  CHARACTER  STCLASS*4,CLHOLD*A,STNU«B*i5,SNHOLD*1 5 

2011  INTEGER  CLASSCMT.TOTCMT.RDCNT.RD 
ion  READ<2O,IOO,EMDs90O)STCLASS,STNUHB,RD 

70H100  F0RHAT(T5,A<,T3,A1S.T116,16) 

BOB  TOTCMTal 

901  CLASSCNT=1 
100«  CLHOLD=STCLASS 

11  OR  SNHOLD>STNUNB 

120R  RDCHT=RD 

130#  ITEHCNT=1 
HO#  URlTE(<i,125) 

1S0R12S  FGRMATCMSTOCK  CLASS  TOTAL  ITEHS  TOTAL  RECURRING  DEHAND") 

UOH  DO  500  K^l, 11000 

170#  READ(20,100,EHD=900)STrLASS,STN0H8,RD 

180#  T0TCNT=T0TCNT+1 

220#  IFISTHUHB.GE.SNHOLD’GO  TO  150 

230#  PRINT,"  ST  NUNBERS  OUT  OF  3E0;REC#  ".TOTCNT 

240#  STOP 

250# 150  SNHOLDsSTNUHB 

240#  IF(STCLASS.EQ.CLH0LD)G0  TO  200 

330#  CLASSCNT»CLASSCNT+1 

340#  URITE<4.)010)CLH0LD,ITEHCNT,RDCNT 

350#  CLHOLD«STCLASS 

340#  RDCNT-0 

370#  ITENCKT*0 

380#200  RDC#T»RDCHT4RD 

390#  IfEHCNT=lTEHCNT+1 

400#500  CONTINUE 

4109900  URITE(4,1010)CLH0LD.ITENCNT,RDCNT 
420#  «RITE(4.950)CLASSCHT 

430#950  FORHATI///"  TOTAL  NUHBER  OF  STOCK  CLASSES  «  ",14) 

440#  «RITE(4,1000)T0TCHT 

450#1000  FORMAT!"  TOTAL  NUMBER  OF  RECORDS  »  ’M4) 

440#1010  F0RHAT(/2X,A4,T18,I9,T33,I11) 

470#  STOP 

480#  END 
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APPENDIX  B 

SELECTION  PROGRAM  FOR  RANKING  BY 
DEMAND  AND  FREQUENCY 
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LIST  GE09 


10  INTE6ER  SC,TD,lARRAY(t31> 

20  CHARACTER41S  SN.HOLOSN 
30  IFLA6=0 

40  H0LDSN=*123454789123454“ 

50  50  REAO(40.100.END=°SOO)SC.SN.IDATE.1DEHAND 
40  100  FORMAtIti, 12, T1.A15,T47, 14, 779,14) 

70  IF(SC.6E.33)G0  TO  200 
90  60  TO  SO 

100  200  1F(SN.NE.HOLDSN)60  TO  300 

110  CALL  TALLY(SN,IDATE,IARRAY,1DEHAND,IT0T) 

120  GO  TO  50 

130  300  1F(1FLAG.E0.0)G0  TO  350 

140  DO  400  K-1,131 

ISO  IBUCKETS=IBUCKETS+IARRAY(K) 

140  400  CORTIMUE 

170  IBUCKETS»IBUCKETS«100000 

180  ITOT»ITOT+f 000000000 

190  UR1TE<S0,425)H0LDSN,IBUCKETS,IT0T 

200  425  F0RHAT(1X,A15,1X,14,1X.I8) 

205  350  IF<SC.GT.53)G0  TO  500 

210  rFLAG=t 

220  H0L0SN=SN 

225  IBUCKETSsO 

230  ITOT=«0 

240  DO  450  J=1.131 

250  IARRAY<J)=0 

240  450  CONTINUE 

270  CALL  TALLY(SN,IDATE,IARRAY,1DENAND.IT0T) 

280  60  TO  50 

320  500  URITE<4,550) 

330  550  FORNATC  END  OF  SN  PROGRAN”) 

340  STOP 
350  END 

340  SUBROUTINE  TALLY(SN.IDATE,IARRAY,IDENAND,ITOT) 
370  CHARACTER  SN«1S 
380  INTEGER  IARRAY<131) 

390  IFdDATE. GE. 4274. AND. IDATE.LE. 9091  )G0  TO  400 
400  URITE(4,S75)SN,IDATE 
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410  575  FORHATC  DATE  OUT  OF  BOUNDS*  ".A15,5X,I4) 

420  RETURN 
430  600  CONTINUE 

440  IFdDATE.GE. 6276. AND. IDATE.LE. 6366) ISUB*(IDATE-626?)/7 
450  IFaDATE.GE.7001 .AND.IDATE.LE.7365)ISUB=(IDATE-6994)/7+13 
460  IFdDATE.GE.OOOl  .AND.IDATE.LE.836S>ISUB=(IDATE-7993)/7-f65 
470  IF(IDATE.GE.9001.AND.IDATE.LE.9091)ISUB=(lDATE-8992)/7+117 
480  lARRAYdSUB)*! 

490  ITOT=ITOT+IDEMAND 
500  RETURN 
510  END 


LIST  S0RT-1B 
0099INN0RN.J 

01004:IDENT:UP1186,AFIT  ACDS  LT  NAIER 

0120I:6HAP:NOECK 

0130!600SH>NOHAP 

0140:SORT:1NOUT,,S 

0tS0!FIELDt(Ct6,C5,C9) 

0160tSEQ:<D2.D3.A1) 

0170:FILCB:1N0UT,*=*,2 

0180:END 

0190$:EXECUTE 

0200$:LINITS:10 

02 1 0$  t  PRHFL : S A , R . S . BOBF /DAT  A 1 -47 

0220«:PRNFL:SZ,U,S.BOBF/DATA1-B 

0230«:FILE:S1,X1R.10R 

0240$:FILE:S2.X2R,10R 

0250«:FILE;S3,X3R,10R 

02604:ENDJ0B 


LIST  SORT-ID 
0099NNNORH,J 

0100t:IDENT;UP1186,AFIT  ACDS  LT  NAIER 

0120«:GHAP:NDECK 

0130!600SH:NONAP 

0140:S0RT:IN0UT,,S 

01S0:FIELD:<C16,CS,C9) 

0160iSE0:(D3,02.At) 

0170:FILCB:1N0UT,»*,2 

OI8O1END 

0190$ {EXECUTE 

0200$:LINITSi10 

0210$!PRNFL;SA,R,S,B0BF/DATA1-47 

0220tsPRHFL:SZ.U,S.B0BF/DATA1-D 

0230$:FILE:S1,X1R.10R 

0240$;F1LE!S2,X2R.10R 

0250$:F1LE:S3.X3R,10R 

0260$:ENDJ0B 
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APPENDIX  C 

PROGRAM  FOR  DATA  ANALYSIS 
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LIST  rtOVElT 


10  CHARACTER413  SN1 ,SN2 

20  INTEGER  BUKET1(130) .BUKET2(130) ,BUKET3(130> 

30  S  REA0(10,10,END=500)SN1,<BUKET1(K>,K-1,130) 

40  READ(20.10,END=i00)SN2,(BUKET2(J),J-1,130> 

30  10  E0RttAT<AtS,5(/23I5)} 

60  IF(SN1.NE.SN2)Ga  TO  700 
70  DO  30  L=1,130 

80  BUKET3(L)>BUKET1<L)4BUKET2(L) 

90  30  CONTINUE 

100  URITE(30.50>SN1,(BUKET3<I),I=1,130) 

110  50  F0RHAT(A15,12I5,7(/15I5)/13I5) 

120  60  TO  5 

130  500  URITE(6,3S0) 

140  550  FORNATC  END  OF  FILE-1") 

150  READ(20,10,END=S75) 

160  URITE<6,590)SN2 

170  590  FORMAT ("  E0F-1,BUT  INFO  STILL  ON  FILE-2, SN2»* ,A15> 
180  STOP 

190  575  URITE(6,595) 

200  595  FORNATC  EOF  ON  BOTH  FILES,  ALL  OK") 

210  STOP 

220  600  URITE(6,660)SN1 

230  660  FORMAT <"  UNEXPECTED  E0F-2,SN1=",A15) 

240  STOP 

244  700  URITE(6,750)SN1,SN2 

246  750  FORNATC  UNMATCHED  SN'S*",A15,2X,A15) 

248  STOP 
250  END 

410  IARRAY(ISUB)»IARRAY(ISUB)-»IDENAND 

420  IF<IARRAY(ISUB).GT.99999)WRITE<6,700)SN 

430  700  FORNAT(  "IDEMAND  HAS  EXCEEDED  99999,  SN=",A15) 

440  RETURN 
450  END 

460  SUBROUTINE  IN1T(IF0UND,IARRAY) 

470  INTEGER  IARRAY(t30) 

480  DO  800  J=1.130 
490  IARRAY(J)>0 
500  800  CONTINUE 
310  IF0UND>0 
320  RETURN 
330  END 
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10  INTEGER  IARRAY<130) 

20  CHARACTER*! S  SN.HOLDSN 
30  CALL  IK1T<1F0UND,IARRAY) 

40  READ<30.10.ENO-500)HOLD3N 
50  10  EORNAT(AIS) 

60  50  REAO<40,100,END-499)SN,10ATE,IOENAND 
70  100  F0RHAT(T1 ,A15,T67,14,T79,I6) 

80  20  IPISN.NE.HOLOSNIGO  TO  130 
90  IFOUND-1 

100  CALL  TALLY(SN,IDATE,IDENAND,IARRAY) 

110  GO  TO  SO 

120  ISO  IFdFOUND.EQ.DGO  TO  300 
130  IF(SN.GT.HOLDSN)GO  TO  200 
140  GO  TO  SO 

150  200  URITE<6.2S0)H0L0SN 

160  2S0  FORNATC*  THIS  3N  NOT  ON  THIS  TAPE=“,A15) 

170  60  TO  450 

180  300  URnE(50,400)HOL0SN.(IARRAY(I),l>1,130) 

190  400  F0RNAT<A15,5(/26I5)) 

200  CALL  INIT(IF0UND,1ARRAY) 

210  430  READaOJO,EHII^SOO)HQLMH 
220  GO  TO  20 

230  499  UR1TE(6,999)H0LDSN 

240  999  FORNATC  END  OF  TAPE=",A15) 

2S0  500  URITE<6,SS0) 

260  350  FORNATC  END  OF  SN  PROGRAH") 

270  STOP 
280  END 

290  SUBROUTINE  TALLY (SN,IOATE,IDEHAND,IARRAY) 

300  CHARACTER  SH«15 
310  INTEGER  IARRAY(130) 

320  IFdDATE.GE. 6276. AND. IDATE.LE.9087)G0  TO  600 
330  URITE(6,S7S)SN,IDATE.IDENAN0 

340  575  FORNATC  DATE  OUT  OF  BOUNDS*  ".AIS.SX.M.SX.IS) 

350  RETURN 
360  600  CONTINUE 

370  IFdDATE.GE. 6276. AND. IDATE.LE.6366)ISUB*dDATE-6269)/7 
380  IFdDATE.GE. 7001. AND. I0ATE.LE.7365)ISUB*dDATE-6994)/7*13 
390  IFdDATE.GE. 8001. AND. IDATE. LE. 8365) ISUB»dDATE-79V3 >77*65 
400  IFdDATE.GE. 9001. AND. IDATE.LE.9089)ISUB»dDATE-8992)/7*117 
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APPENDIX  D 

PROGRAM  TO  COMPUTE  FORECASTS  AND  PERFORM 
COMPARATIVE  ANALYSIS 
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lll^C 

lZf=C 

13fsC*«i 

STAR  PROGRAfl  CORPUTES  FORECASTS  USIW  IHJLTIPLE  TECHAIBUES 
AND  PERFORRS  BASIC  CORPARATIVE  ANALTSIS 

14l> 

PROCCAR  STAR(INPUT.OUTPUTiTAPE10*INPUT,TAPE6*OUTPUT) 

15f: 

CORRON  NSN (Z) iNDAT ( 150) i IFSE(1S0) . IFDE(1S0I • 

1M= 

1IFAE(1S0) .IFRA4(1S0) i IFRA8(1S0) f 

171* 

lRS(130)>INR(150liISEQ(3500l 

181* 

REAO(10>9ai)  I 

191*901 

FORRATIU) 

ZIf* 

L*I 

Zll* 

DO  899  J*1,Z5 

ZZI* 

READ ( 101 900) ISEQ ( I ) iNSN ( 1 ) .NSN (Z) . 

Z3f* 

l(NDAT(IZ).IZ*lilZ).lSEa(m).  INDAT 

Z4|* 

1 (IZ) .IZ*13iZ7) . ISE8 (l^Z) f (NDATIIZ) . 

Z50* 

lIZ*Z8.4Z).ISE0(U3)i(NBATIIZIf 

Z6I* 

lIZ*43.57).lSE0(M)t  (NDATIIZ). 

Z70* 

1IZ*58.7Z).ISEQ(I>5). (NDATIIZ). 

Z80* 

1IZ*73.87)  .ISEOIIH) .  (NDATIIZ)  > 

Z9#= 

lIZ*88il0Z).ISEQ(I>7). (NDATIIZ) I 

300* 

lIZ*103.U7).ISEQ(l4a). (NDATIIZ). 

310* 

117*118.130) 

3Z0*900 

FORRATIU, 1I.A10.A5,1ZI5.7I/I4. 1515). /I4.13I5) 

330* 

J5*L+8 

340* 

DO  105  J1*L.'J5 

350* 

IFIISEQIJD.NE.JI)  CO  TO  898 

360*105 

CONTINUE 

370* 

CALL  SINCEI 

380* 

CALL  OOUBEI 

390* 

CALL  ADAPEX 

400* 

CALL  NOVA 

410* 

CALL  ROV8 

4Z0> 

CALL  RULTI 

430* 

CALL  ERROR 

440* 

L*L*9 

450* 

i*L 

460*899 

CONTINUE 

470* 

STOP 

480*898 

PRINT  975,ISEQIJ1-1) 

490*975 

FORRAT(5X,Z0KCARD  OUT  OF  SEOUENCE.AI.U) 

5N* 

STOP 

510* 

END 
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53f*C  SUBROUTINE  SINCEI  CONFUTES  SINGLE  EIFONENTIAL  FORECAST 
S4f«C  VALUES 

S6f«  SUBROUTINE  SINGE! 

57B-  CONNON  NSN(2)iNDAT(lSg)iIFSE(lSfliIFI)E(15f)> 

58f*  lIFAE(lSB)iIFNA4(lSf)iIFNA8(lM>> 

59|:  INS ( 131) I  INN ( ISf ) I ISEQ (35(B) • 

IFSE(l)-NOAT(l) 

&tf>  00  m  UMZ9 

iZB-  FSEMFSE(I)M.Z*(NDATIII-IFSE(I))) 

63B=  IFSE(I+1)*IFI!(FSE) 

&4B>Z(B  CONTINUE 
65B:  RETURN 

MB-  END 

&7(-C 

49(<C  SUBROUTINE  DOUSE!  CONFUTES  DOUBLE  E!FONENTIAL  FORECAST 

7BB*C  VALUES 

71B*C 

72B:C**«M«*******M***«M«****4**M*tM«4t«*t*tt4lt(»««4t»ff««t»tf 

73|:  SUBROUTINE  OOUBE! 

74B«  CONNON  NSN(Z)fNDAT(15B))IFSE(lSB)fIFDE(lS()> 

75B:  1 IFAE ( tSB) > IFNA4( ISB) f IFNA8 ( 15B) • 

7&B>  lNS(13B))INN(tSB)!lSEa(35BB) 

77B*  IFDE(1)»NDAT(1) 

7«#» 

79B^  STC»B.B 

8BB>  DO  3BB  I >2 i 129 

81B*  ATC*1FSE(I)-IFSE(I-1) 

82B«  STC=STC+(.24(ATC*STC)) 

S3B»  IF(STC.LT.B.)  STC*B.B 

84B»C  5«l/.2 

85B«  FDE»IFSE(I)4(54STC) 

8«B«  IF0E(I41)>IFI!(FDE) 

87(<3BB  CONTINUE 
88B<  RETURN 

89B>  END 
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9ZfcC  SUBROUTINE  MMPEI  CONPUTES  ADAPTIVE  EIPONENTIAL  FORECAST 

93f:C  VALUES 

94I*C 

9il=  SUBROUTINE  AOAPEI 

97l<  COmON  NSN(Z)iNDAT(lSI))IFSE(lSII.IFDE(15f). 

98f:  lIFAE(lSI)>IFNA4(lSf)>IFnA8llS0l. 

99f-  lNS(13lliinn(lSf)>ISEQ(3Sff) 

IBM:  IFAE(ll<NOAT(U 

nil::  ERTtl.f 

1IZB>  EKRT^.f 

1I3B°  00  411  I-lflZ9 

1B4B=  IET^NDAT(I)-IFA£(I) 

1150^  IFdET.EO.I.I  CO  TO  371 

1060:  CO  TO  380 

1070:370  ALPHA:0.i 

1080:  ERT:0.0 

1090:  EHRT:0.0 

1100:  CO  TO  390 

1110:380  ERT:(.Z*IET)+(.8*ERT) 

1120:  IET:IABS(IET1 

1130:  ENRT:(.2»IET)  +  (.8»EIWT) 

1140:  ALPHA:AB3IERT/EHfiTl 

1150:390  FAE:IFAE ( I) ♦ (ALPHA* (NDAT ( I ) - IFAE 1 1) ) » 

1140:  1FAE(I*1):IFIX(FAEI 

1170:400  CONTINUE 

1180:  return 

1190:  end 

1Z00:C 

1Z10:C************************************************** ************ 

1ZZ0:C  SUBROUTINE  Ni}V4  CONFUTES  4  TERN  NOVINC  AVERACE 

1Z30:C  FORECAST  VALUES 

1Z40:C 

1Z50:C************************************************************** 

1Z60:  SUBROUTINE  N0V4 

1Z70:  CONNON  NSN(2)>NDAT(150)rIFSE(150)rIFDE(150>r 

1Z80:  IIFAE ( 150) > IFNA4 (150) i IFNA8 (150) > 

1Z90:  lNS(130)fINN(150)iISEQ(35M) 

1300:  DO  490  Il:ll4 

1310:490  IFNA4(II):N0AT(I1) 

13Z0:  IF)M4(1):N3AT(4) 

1330:  DO  SM  I:4tlZ9 

1340:  niA4: (NDAT ( I ) ♦NOAT ( I- 1 ) tNOAT (I-Z) *NDAT (1-3) ) 74 

1350:  IFNA4(I*1):IFIX(FNA4) 

1360:500  CONTINUE 

1370:  return 

1380:  end 
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141f<  SUBROUTINE  IIOV8  COdPUTES  8  TERN  ROVING  AVERAGE 

UZI^C  FORECAST  VALUES 

1A3I< 

1458=  SUBROUTINE  N0V8 

14&f:  CONRON  NSR(2)>NDAT(lSII>IFSE(15f)iIFDE(lSI)i 

1478^  lIFAE(lS8)iIFRA4llS8)>IFRA8(158li 

1488^  lHS(138)iIRRI158).ISE0(3S88l 

1498==  DO  598  I1M)8 

1588^598  IFRA8(I1)=NDAT(I1) 

1518=  IFRA8(1)=NDAT(8) 

1528=  00  488  I =8 1129 

1538=  FRA8=(NDAT(I|4N0AT(I-1)+NDAT(I-2I4NDAT(I-3I+ 

1548=  1N0AT(I-4)4NDAT(I-5)4N0AT(I-&)4NDAT(I*7))/8 

1558=  lFRA8a+l)=IFIX(FRA8> 

1548=688  CONTINUE 

1578=  RETURN 

1588=  END 

1598=C44444444444444444444444444444444444444444444444444444444444444 

1488=C  SUBROUTINE  IHILTI  CWPUTES  RULTIPLE  RODEL 

1418=C  FORECAST  VALUES 

1428=C 

1438=C44444444444444444444444444444444444444444444444444444444444444 

1648=  SUBROUTINE  NULTI 

1458=  CORRON  NSN(2)iN0AT(l58)fIFSE(lS8)tIFDE(158lf 

1468=  UFAEU5I1  .IF«ft4U58WIFRft8U5«> . 

1478=  lHS(138)iIRR(lS8lfISEai3588> 

1488=  DO  498  Il=lf9 

1498=498  NS  (ID  =8 

1788=  DO  728  I=9il29 

1718=  NVAR8=N0AT ( I » *8087 ( I- 1 ) +8087 ( I-2)+NDAT(I-3) 

1728=  NVAR1=(NVAR8-(IFSE(1)*IFSE(I-1>+1FSE(I'2)+ 

1738=  lIFS£(I-3))) 

1748=  NVAR1=IABS(NVARD 

1758*  NVARZ=IABS(NVAR8-(lF0E(D4lFDEII-U4lFDE(I-2)4 

1748*  lIFDE(I-3))) 

1778=  NVAR3=IABS(NVAR8-(IFAE(D4lFK(M)4lFI€tI-2)4 

1788*  lIFAE(I-3))) 

1798=  NVAR4=IABS(NVAR8-{IFRA4(II4IFRA4(I-1)«!FNA4(I-2K 

18H=  lIFNA4(I-3))) 

1818*  NVAR5* I ABS (NVAR8- ( IFRA8 ( I ) 4IFRA8 (1-1)4 

1828*  lIFHAe(I-Z)4lFNA8(I-3))l 
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B3I* 

M*1 

184f* 

NVAR*NVAR1 

185f* 

IF(NVARZ.LT.NVAR)  CO  TO  785 

1841*711 

IF(NVAR3.LT.NVAR)  GO  TO  784 

187I*7BZ 

IF(NVAR4.LT.NVAR)  CO  TO  787 

188«*7I3 

IF (NVAR5.lt. NVAR)  CO  TO  788 

1891= 

CO  TO  718 

19M*7i5 

H*Z 

1911* 

NVAR*NVAR2 

1921= 

CO  TO  781 

193f*7f4 

11*3 

1941* 

NVAR*NVAR3 

1951* 

CO  TO  782 

194i*7r 

M*4 

1971= 

NVAR*NVAR4 

1988* 

CO  TO  783 

199f*7f8 

H*5 

2««B*71B 

MS(I)*H 

2818* 

IF(NS(I).EQ.l) 

IHM(Id)*IFSE(l4l) 

2828* 

IFinSdl.EQ.Z) 

INHddtsIFDEdd) 

2838* 

IF(HS(I).EQ.3) 

IHM(l+l»®IFAEdd» 

2848* 

IF(RS(I).EQ.4) 

1HM(I41»*IFHA4(I+11 

2858* 

IF(HS(II.EQ.5) 

INR(l41)*IFMA8d4U 

2848*728 

CONTINUE 

Z878* 

RETURN 

2888* 

END 

Z1M<  SUBROUTINE  ERROR  MEASURES  THE  DECREE  OF  ERROR  AND  BIAS  IN 

ZUf*C  THE  FORECASTS 

ZUf<C444H444444444444444444444444444444 444444 44 444444444444444444444 

Zl3f>  SUBROUTINE  ERROR 

Zl*9>  COmON  NSNIZIfNDATdSBIflFSEilSBIiIFDEdSBW 

Z1 IIFAE  dSf  I .  IFRA4  (151)  >  IFHA8  dSB)  i 
Z16f=  lHSd3l)dmd59)tISEQ(35M) 

Zl7f‘  DIMENSION  ICTRdSf) 

ZlBf:  SEHAD<f.l 

Z19B*  DEnAD<l.f 

ZZ«f:  AEHAD>I.I 

ZZlf'  A4MAO*l.f 

ZZZI>  A8MA0*f.f 

ZZ3f>  SERE*f.f 

Z24I*  DERE*l.f 

ZZ5I>  AERE<l.f 

ZZ4B<  A4RE<B.f 

ZZ7f>  A8RE*f.f 

ZZ8f>  SEPE*f.l 

Z2fB<  OEPE*l.l 

23M>  AEPE*l.f 

231f>  A4PE*B.f 

2321*  A8PE>IJ 

233B*  FMB*I.I 

234f*  FRE*f.f 

235#*  FPE*I.I 
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36#»  DO  8M  1*13.131 

237#*  ERR  s  RDAKD-lFSEd) 

238#*  SERAD*#8S(ERR»+S£I(«) 

239#*  SERE*ERR+SER£ 

29##*  IFCERR.EQ.I.)  CO  TO  757 
2«1#*  IFCRMKD.EO.I.)  CO  TO  731 
2«2»*  SEPE*  (ABS  (ERR/NOAT  ( 1 )  I  •!«.  »+SEPE 

243#=  CO  TO  752 

244#*75I  SEPE*(ABS(ERR/1)#1I#.»+SEPE 

24S#*752  ERR*NBAT(1)-IF0£(1) 

244#*  OEMD*ABS(ERR)40ERAO 

247#*  DERE*ERR+OERE 

248#=  IFCERR.EO.i.)  CO  TO  754 

249#*  IF(iniAT(l).EQ.|.l  CO  TO  733 

25##*  OEPE* (ABS (ERR/BOAT ( 1 ) I «tf|. ) tOEPE 

251#*  CO  TO  754 

252#* 753  DEPE*  (ABS  (ERR/1 )  *11# .  HDEPE 

253#*754  ERR*l(DAT(l)-IFA£(n 

254#*  A£»(AD*ABS(ERR)*AE«AD 

255#*  AERE=ERR*AERE 

254#*  1F(ERR.E0.#.|  CO  TO  754 

257#=  IF(NDAT(I).EQ.i.)  CO  TO  755 

258#*  AEPE* (ABS (ERR/MDAT ( I ) ) *1##. )+A£PE 

259#*  CO  TO  754 

24##*755  AEPE*(A8S(ERR/1I*1|#.)+aePE 

241#*754  ESR*t(0AT(l)-IFHA4(I) 

2421*  A4«AD*ABS(ERR14AWAD 
2431*  A4R£.ERR<A48£ 

244#*  lF(ERR.EQ.f.l  CO  TO  758 
245#*  lF(l(DAT(l).Efl.i.l  CO  TO  757 

244#*  A4PE*(A8S(ERR/MDAT(I))«1#|.)+A4PE 

247#=  CO  TO  758 

Z48#*757  A4PE*(ABS(ERR/1I»1I#.)+#4PE 

249#*758  ERR*raRT(I)-IFHA8(I) 

27##*  A8fMD*ABS(ERR)*A8nA0 
271#*  A8RE*ERR*A8RE 

272#*  1F(ERR.EQ.#.1  CO  TO  74# 

273#*  IF(W1AT(1).EQ.#.)  GO  TO  739 

274#.  A8PE*(ABS(ERR/RDAT(in»l##.)*A8PE 

275#*  CO  TO  74# 

274#*759  A8PE*(ABS(ERR/1)*1##,)*A8PE 

277#*74#  ERR«IIDAT(1)-1IW(1) 

278#*  FIIAD*ABS(ERR)*F)MO 

279#*  FRE*ERR*fRE 

28##*  1F(ERR.E8.#.»  CO  TO  745 

281#*  lF(iaiAT(l).E8.#l  CO  TO  741 

282#*  FPE*(ABS(ERR/W)AT(l))ttN.)*FPE 

283#*  CO  TO  745 

284#*741  FPE*IABS(ERR/1)*1H.)+FPE 

28S#*745  continue 

284#*###  CONTINUE 
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•7i*  SEBA0*SEIWD/U8. 

SCRE*(SERE/U8.>»(-1.») 

28«*  SEPE«SEPE/118. 

MM*  0ENA0*0EIW0/U8. 

Mll=  DERE*(DERE/lI8.)«(-l.f) 

M2«*  DEPE»DEPE/118. 

M3I*  AEHA0*AEM0/U8. 

M4I=  AERE=(AERE/118.)»(-1.I» 

M5I*  AEPE*AEPE/ll8. 

M6I=  A4RAD*A4M0/118. 

M7I=  A4R£s1A4RE/118.)»(-1.II 

M8I=  A4PE*A4PE/118. 

M9I=  A8HAD*A8nAD/118. 

38M*  A8fi£=(A8fi£/U8.)«(-lJ| 

Mil*  A8PE=A8PE7118. 

3KI*  FHA0*FHAB/n8. 

3131=  FfiE»(FRE/118.)»(-l.g) 

3841=  FPE=FPE/1I8. 

385#*  ICTR(1»*1 

3841=  00  848  1=2,138 

3878=848  ICTR(I)=ICTR(I-1I+1 

3888=  l«ITE(4,9a8)  NSNd) ,NSN(2I 

3898=988  F0RRAT(1H1,///51,5HNS»;  ,2A18) 

3188=  13*1 

3118=  14*15 

3128*  II«i 

3138*  00  878  12*1,8 

3148*  IF(II.£0.4>  CO  TO  845 

3158*  CO  TO  841 

3148*845  HRITE(4,912) 

3178=912  FORHATdHl) 

318#*  11*1 

3198=841  l«nE{4, 9811  (ICTR(I), 1=13, 141 
3288=981  FORRAT  (//5I  ,8H»»«£EKm,  41 , 13, 14I7J 

3218=  URITE«4,9821  (MOATU)  ,1*13,141 

3228*982  FOfinAT(/53,4HACTUAL,4J, 15, 14171 
3238=  WITE(4,9g3l(IFSE(I},  1=13, 141 
3248=983  F0R»1AT(/5I,8H1-SIWLE,2M5,14I71 
3258=  »«1ITE{6,984I(IFDE{I1, 1=13,141 

3248*984  FORRAT (/5I ,8H2*D0UBLE,2l, 15, 14I7J 
3278=  IIR1TE(4,9851(IFAE(I), 1=13,141 

3288*985  F0RIWTt/5J,7H3-ADAP7,3I,I5,l4I7J 
3M8*  l«ITE(4i9841(IFHM(n,I*l3,I4) 

3388*984  F0R«AT(/5J,7H4-4TER«,3M5, 14171 
3318*  HRITEIO, 9871  lIFlVWd), 1*13,14) 

3328*987  F0R|IAT(/5I,7H5-8TERR,3I,I5,14I71 
3338*  WIITE(4, 9881 (HSll), 1*13, 141 

3348*988  F0R1IAT(/5I,4HBETH0D, 41,15,1417) 

3358*  WITE(4,989>(I)))II1),I*I3,14) 

3348*989  FOWttT(/5l,8H)«JLTIBOO,2I,I5,14I7) 

3378*  13*13+15 

3388*  14*14+15 

3398*  II.IUi 

3488*878  CONTINUE 
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411*  URlTE(4f9fl)(ICTR(I)<I*I3il3l) 

3421*  WITE(4>91f) 

3431*911  F0WflT(lH-»i9IIf3HnADt8Ii4HBlAS«Mi4HMP£) 
3441*  URITE(4i9flZI(NMT(IltI*[3>13f) 

3451*  l«ITE(4.9g3)(IFS£U)il*I3»13l) 

344f*  URITE(4>911)  SENADiSEREiSEPE 

3471*911  F0RnAT(lH4.84Xi3F11.4) 

34M*  URlTE(6i9fl4)(IF0E(I)iUI3fl3l) 

3491*  URITE(4i911)  DEHAOiOEREiOEPE 

3SM*  URITE(4i90S)aFA£(niI*I3fl3f> 

3511*  1»R1TE(4»911)  AEIlAD.A£REiA£PE 

3521*  IIRITE(6.9«41  (IFHA4(I)  it*I3il3«) 

353**  »RITE(4.911)  A4l(ADiA4RE.A4P£‘ 

3541*  BRlTE(4i9i7)(IFRA8mit=I3il3«» 

3551*  IIRITE<4i911»  A8«A0fA8R£fA8P£ 

3541*  IIRITE(4>9R8>(IIS(n>l'I3il3l) 

3571*  MRITE(4t989KIWI(I)>l=I3.13l» 

3588*  8RITE(4.9111  FRAO.FREiFPE 

3591*  R£TURI» 

34H*  £NO 
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